Localisation for non-monotone Schrodinger 

operators 

Alexander Elgart 1 , Mira Shamis 2 , and Sasha Sodin 3 
October 16, 2012 



43 

Abstract 

We study localisation effects of strong disorder on the spectral and dy- 
namical properties of (matrix and scalar) Schrodinger operators with non- 
monotone random potentials, on the ci-dimensional lattice. Our results in- 
clude dynamical localisation, i.e. exponentially decaying bounds on the tran- 
sition amplitude in the mean. They are derived through the study of frac- 
tional moments of the resolvent, which are finite due to resonance-diffusing 
effects of the disorder. One of the byproducts of the analysis is a nearly opti- 
mal Wegner estimate. A particular example of the class of systems covered 
by our results is the discrete alloy-type Anderson model. 



1 Introduction 

1.1 Random Schrodinger operators 

The prototypical model for the study of localisation properties of quantum states 
of single electrons in disordered solids is the Anderson Hamiltonian Ha, which 
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acts on £ 2 (z d ) by 



m adjacent to n 

where the entries u(n) of the potential are random and independent. 

The basic phenomenon, named Anderson localisation after the physicist P. W. 
Anderson, is that disorder can cause localisation of electron states, which mani- 
fests itself in time evolution (non- spreading of wave packets), (vanishing of) con- 
ductivity in response to electric field, Hall currents in the presence of both mag- 
netic and electric field, and statistics of the spacing between nearby energy levels. 
The first property implies spectral localisation, i.e. the spectral measure of Ha is 
almost surely pure point, and almost sure exponential decay of eigenf unctions. 

These properties are known to hold for Ha in each of the following cases: 1) 
high disorder (the coupling constant g is large), 2) extreme energies, 3) weak dis- 
order away from the spectrum of the unperturbed operator, and 4) one dimension, 
d= 1. 

Historically, the first proof of spectral localisation was given by Goldsheid, 
Molchanov and Pastur [fT2l . for a one-dimensional continuous random Schrodinger 
operator. 

In higher dimension, the absence of diffusion was first established in 1983 by 
Frohlich and Spencer using multi-scale analysis. Their approach has led to a 
multitude of results on localisation for a wide range of problems. The reader is re- 
ferred to the monograph of Stollmann ifTTll or the recent lecture notes of W. Kirsch 
ffT4l for a review of the history of the subject and a gentle introduction to the multi- 
scale analysis — which is not used here. 

One of the ingredients of multi-scale analysis is the regularity of the integrated 
density of states, the (distribution function of the) average of the spectral measure 
over the randomness. 

Ten years later Aizenman and Molchanov [2] introduced an alternative method 
for the proof of localisation, known as the fractional moment method, which has 
also found numerous applications. In particular, in UJ, Aizenman introduced the 
notion of eigenfunction correlator, which, combined with the fractional moment 
method, allowed him to give the first proof of dynamical localisation. We refer 
to the lecture notes of Stolz |[T9ll and Aizenman and Warzel [3] for a survey of 
subsequent developments. 

In the fractional moment method, an a priori estimate on the diagonal elements 
of the resolvent (Ha — A) -1 plays a key role in the underlying analysis. 



2 



In many situations, regularity of the integrated density of states follows from 
the regularity of the distribution of the potential. This was first proved by Wegner 
ll20ll . therefore regularity estimates on the density of states are called Wegner esti- 
mates. An essential ingredient in his argument is the monotone dependence of the 
spectrum of Ha on the random variables v(n). A modification of this argument 
was applied by Aizenman and Molchanov to give an a priori bound on the average 



The monotone dependence of the spectrum of the random variables is also 
used in the fractional moment's proof of dynamical localisation (via a variant of 
spectral averaging). 

Recently, several challenging problems (listed below) arose in different con- 
texts, in which the dependence of the spectrum on the random variables is not 
monotone. In this paper, we develop a strategy to prove localisation which is 
applicable to some of these models. 

One close relative of the original Anderson Hamiltonian H A is the random alloy- 
type model, in which the potential V(n) at a site n 6 Z d is obtained from inde- 
pendent random variables i>(m) via the formula 



where the index k takes values in some sub-lattice T of Z d . If all the coefficients 
ak have the same sign, the system is monotone, since the dependence of the spec- 
trum on V is monotone. Localisation in such systems is well understood by now, 
even in the continuum setting. The existing technology is however not well suited 
to the non-monotone case, i.e. when a k are not all of the same sign. Mathemati- 
cally, the problem becomes especially acute when ^ a k = 0. 

There is no physically compelling reason for a random tight binding alloy 
model to be monotone, and the natural question is whether Anderson localisation 
still holds if one breaks the monotony. 

Non-monotone models also naturally appear in the class of block operators. In one 
such model, introduced by Frohlich, and studied by Bourgain in [4], the matrix- 
valued potential is given by V(n) = U (n)*AU (n), where A is a fixed self-adjoint 
r x r matrix, and U(n) are independently chosen according to the Haar measure 
on SU(r). Bourgain proved a volume-dependent Wegner estimate and Anderson 
localisation near the edges of the spectrum using methods from complex analysis. 



of \{H A - A-iO)- 1 ^)! 5 . 




(1.1) 
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Another class of non-monotone random block operators, associated with the 
Bogoliubov - de Gennes symmetry classes, was studied by Kirsch, Metzger, and 
Miiller in [fT5l and by Gebert and Miiller in ifTOl . The former paper establishes the 
Wegner estimate for such random operators in a neighborhood of the inner band 
edges, by recovering some form of the monotonicity. The latter paper uses the 
bootstrap multi-scale analysis to prove dynamical localization in the same regime. 

The original motivation for this work was to study a problem suggested by Tom 
Spencer, in which the matrix- valued potential is of the form 

V(n) = ( ^ %), 
y ' \ a — v(n) J 

where the variables v (n) are independent and identically distributed. If the distri- 
bution of v(n) has bounded density, the eigenvalue distribution of a single V(n) 
is 1/2-Holder; this is optimal, since the density of the eigenvalue distribution di- 
verges as | A =F a\~ l l 2 at the energies ±a. Spencer conjectured that the integrated 
density of states for the full Hamiltonian is also at least 1/2-Holder. 

Electromagnetic Schrodinger operators with random magnetic field, the random 
displacement model, random block operators and Laplace-Beltrami operators with 
random metrics are other examples of systems with non-monotone parameter de- 
pendence which were intensively studied recently. 

We refer to the paper of Elgart, Kriiger, Tautenhahn, and Veselic for a 
survey of recent results on the sign indefinite alloy-type models with T = Z d 
and the bibliography pertaining to some of the models mentioned in the previous 
paragraph. 

Summary of results. For the last few years there has been a continuous effort 
to bring the understanding of models with non-monotone dependence on the ran- 
domness to a same level as the one for monotone models. In this paper, we present 
a method to prove Anderson localisation and a Wegner estimate for several non- 
monotone models, achieving this goal. Theorem ll.il and its corollaries pertain to 
a class of models with matrix- valued potentials; when applied to Spencer's model, 
it shows that, under some assumptions on the distribution of v (n), the integrated 
density of states is (1/2 — e)-H6lder for any e > 0, at large coupling g. Unfortu- 
nately, Theorem [Tj] does not directly apply to Frohlich's model. 

Theorem 1 1.21 and its corollaries establish Anderson localisation and a Wegner 
estimate for the alloy-type model (ll.lt , in the case that a is finitely supported. 
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Our argument can be viewed as a further augmentation of the fractional mo- 
ment method of Aizenman-Molchanov [2J. In particular, Proposition 12.11 is a 
modification of [2, (2.25)], whereas Proposition B.ll is a version of the decoupling 
estimates [2, Lemmata 2.3,3.1]. The innovation of this work is that we do not rely 
on an a priori estimate on the moments of diagonal resolvent elements; instead, we 
prove such an estimate in parallel with localisation. We also propose an argument 
which allows to deduce dynamical localisation directly from the resolvent esti- 
mates, and which works in the non-monotone setting as well as in the monotone 
oneQ 

Relation to some past and present works. For the one dimensional continuum 
alloy-type random models the proof of the complete Anderson localization was 
first given by Stolz Ifl8l 

Outside the spectrum of the unperturbed operator (corresponding to the ran- 
dom potential being switched off) one can obtain Lipschitz regularity of the inte- 
grated density of states by reducing the problem to the monotone case. The opti- 
mal Wegner estimate in this case was established by Combes, Hislop, and Klopp 
AH (in the continuum, but their argument is equally applicable in the discrete set- 
ting). This input can be used to prove Anderson localisation in the regimes of 
extreme energies and weak disorder away from the spectrum of the unperturbed 
operator (for the latter regime for the continuum models this result goes back to 
Klopp OH). 

Our work covers the remaining perturbative setting, namely the high disorder 
regime, where we prove complete localisation. 

Recently, Bourgain (private communication) devised a different approach that 
allows to prove s-regularity of the density of states for a wide class of non- 
monotone models which includes Frohlich's model, as well as some of the models 
we consider in this note. 

1.2 Notation and statement of results 

Let Q = (V, 8) be a graph with degree at most k; the set of vertices (sites) V 
may be either finite or countable. The main example is the lattice Q = Z d (where 
k = 2d), however, the greater generality does not require additional effort here. 

'After this work was completed, we have learned from Victor Chulaevsky that a similar argu- 
ment was used earlier by Germinet and Klein in the context of the multi-scale analysis, see ifTTl 
Theorem 1 , Corollary 1 ] . 
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For x, y G V, denote by dist(x, y) the length of the shortest path connecting x to 
y; when Q = Z d , 

dist(z,y) = \\x - 2/|| i . 

Let 

w : Q x V — >■ R 

be a collection of independent, identically distributed random variables, where 
(fi, P) is a probability space, and we assume that the distribution /i of every v(x) 

Al is a-regular for some a > 0, meaning that n\t — e, t + e] < Caic" for any 
e > and t G R; 

A2 has a finite g-moment for some q > 0, meaning that / < Ca2- 

For example, the Gaussian distribution and the uniform distribution on a finite 
interval satisfy Al with a = 1 and A2 with any q > 0. 

We shall denote the expectation by (•) and the expectation over the distribution 
of onev(x) by (-) v ( x ). 

In the electron gas approximation the system of electrons in a crystal is mod- 
eled by a gas of Fermions moving on a lattice. The excitations of the system are 
described by an effective one-body Hamiltonian H, which consists of a short- 
range hopping term and a local (single site) potential. Each site x of the lattice 
will be assumed to have k internal degrees of freedom. 

Single site (matrix) potential: For any x G V, define a Hermitian matrix 

V(x) = v(x)A(x) + B(x) , 
where the Hermitian k x k matrices A(x) and B(x) satisfy 

bi po^iiJAOr)- 1 !! <c B i; 

B2 \\B{x)\\ < CW 

Hopping: For every ordered pair (x, y) G V x V of adjacent sites (i.e. (x, y) G £) 
we introduce a /c x k matrix (kernel) K(x, y) so that 

B3 K(y,x) = K(x,y)* and y)|| < C B3 . 
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We are now in position to introduce our one-particle Hamiltonian. Namely, 
let if be a random operator acting on £ 2 {V) <8> C k (the space of square-summable 
functions ip : V — > C fc ) 

(Hil>)(x) = V(x)<4>(x) +g- 1 Y,K(x,yW(v) , (1-2) 

where g > is a coupling constant, and the sum is over all y E V such that 

(x,y) E S. 

Let G\ = (H — A) -1 be the resolvent of H, A ^ R. It is known that the 
limit G\ + io = lim e _^ +0 G\ + i e exists for almost every A 6 R. In the following, 
(\\G\ +i0 (x, y) can a priori be formally interpreted as 

Km(\\G x+ie (x,y)\\ s )-, 

e— >-+u 

a posteriori, G\ +i0 is finite almost surely, and 

lim (\\G x+ie (x,y)\\ s ) = (II lim G x+i€ {x,y)\\ s ) . 

Theorem 1.1. Let < s < 2 ka+kq ' There exists C > f/za? mav depend on a, q, 
Cai-Cb3 s such that for any A G R and any g > Ck 1 / 3 / (1 + |A|) 

C ( Ck \ dist ( x ' y ^ 

(\\Gm*(x,v)V) < { g s {1 + l x\y ) 

Let us state some corollaries for the homogeneous setting; for simplicity, assume 
that Q = Z d (we denote the vertices of Z d by m, n, • • • ) We also assume that 

C A(m) = A, .B(m) = B, K(m, n) = K(m - n). 

The density of states p is defined as the average of the spectral measure corre- 
sponding to H: 

J /(A)dp(A) = itr</(£0(n,n)>, 

where tr stands for the trace. The integrated density of states is the distribution 
function A M- p{— oo, A] of p. 

The assumption C guarantees that these definitions do not depend on the 
choice of the vertex n £ Z d . 

Theorem II . 1 1 implies the following Wegner-type estimate: 
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Corollary 1.1.1. Assume C. If g > Cd 1 ^ / (1 + |A|), then the integrated density of 
states is locally s-Holder at Xfor 

aq a 
2ka + kq k (l + — j 

uniformly in g -)■ oo: 

p[A-e,A + e] < C(l + |A|)- S e s . 

In particular, for any distribution with bounded density and finite moments the 
integrated density of states is l/(k + e)-H6lder for any e > 0. 

Next, Theorem 11.11 implies dynamical and spectral Anderson localisation: 

Corollary 1.1.2. Assume C. Let I be a finite interval of energies, and let 

Cd 1 ' 3 

9 > • 

1 + min Ae/ |A| 

Then, for any m ^ n G Z d , 

s dist(m,n) 

(sup|e^(m,n)|) < Cdist(m ; n) M ( a( , C * a ) ^ , (1.3) 

where Hj = PjHP^ Pj is the spectral projector corresponding to I. Therefore 
the spectrum of H in I is almost surely pure point. 

The first part of the last corollary follows from Theorem ll.li (14.11) . and Theo- 
rem [4T2J The "therefore" part follows from the summability of the right-hand side 
of (11.31 ) via the Kunz-Souillard theorem (6l Theorem 9.21]. 

1.3 Extensions 

Alloy-type models. Consider the operator H with potential (11.11) acting on t 2 ( Z d ) . 
Let B n be the set of v(m) for which a n _ m 7^ 0. We will impose the following 
assumptions on the random potential: 

1 . the set B n is non empty for all n; 

2. the cardinality k = #{m | a m ^ 0} < 00; 
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3. the distribution of v(ta) satisfies Al and A2. 



Theorem 1.2. Let < s < -, 

and any g > Cd 1/s /(l+ |A|) 



aq 



. There exists C > such that for any A G R 



2ka+kq 



C ( Cd 



) 



dist(m,n) 



(|G A+J0 (m,n)| s )< 



(1 + |A|)« U S (1 + |A|) 



Similarly to Corollaries I1.1.1I1.1.21 one can deduce a Wegner estimate and 
Anderson localisation. 

Relaxing the covering condition. The assumption Bl is usually referred to as 
a covering condition. In our analysis, it enters in the proof of Lemma IZ2l In 
particular, all our results are still valid (albeit with the less sharp bound on the 
underlying localisation length) if one replaces Bl with 



For a fixed non zero matrix A(y) and a generic matrix B(y) the estimate Bl' holds 
true for g large enough. For instance Bl' is applicable (with a = 0) for 
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in particular, for suggesting to use eigenfunction correlators in the proof of dy- 
namical localisation, to Giinter Stolz for remarks on a preliminary version of this 
paper, and to Victor Chulaevsky for bringing the work of Germinet and Klein 0TJ 
to our attention. 

2 Proof of theorems 

In the proof of Theorem ll.il we assume that the graph Q is finite and A ^ R. The 
estimates will be uniform in #V — > oo denotes cardinality) and Im A — > 0, 



BV (\\(V(y) - Xy'W*}^ < Cg a ,witha < s. 
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therefore the statement for infinite graphs and real A can be deduced as follows. 
First, an infinite graph can be approximated by its finite pieces; the matrix ele- 
ments of the resolvent corresponding to the finite pieces converge to the matrix 
elements of the resolvent corresponding to the infinite graph, yielding the same 
estimate for A ^ R. Then one can let Im A go to zero. 

The proof of Lemma 12 .21 below will be postponed until Section[3j 

Proposition 2.1. For any s < 2 ka+kq there exists C > (depending on s and the 
constants in the assumptions ) such that for any A ^ R 

C 

2(1 + |A|) 



Gx(x,y)\\ s ) < f { g- s J2(\\Gx(x,z)\n+5 : 



xy 



where 



&xy 



1, x = y 
0, x^y 



is the Kronecker 5. 

Proof. By definition of G x , 

G x (x, y)(V(y) - A) = -g' 1 £ G x (x, z)K(z, y) + 5 xy . 

Therefore 

(\\G x (x,y)(V(y) - A)|| s ) < <?«,(,- ^(||G A (x, z)\\ s ) + 5 xy . 

Lemma 2.2. For s < 2 ka+kq ' there exists C ( depending on s and the constants in 
the assumptions) such that 

(\\G x (x,y)(V(y) - A)|| s ) > C'^WG^x, y)\\ s ) (1 + |A|) S . 

The proposition follows. □ 

Corollary 2.2.1. For any s < 2fc °^ fc > we have 

max(\\G x (x,y)\\ s ) = (\\G x (x,x)\\ s ) , 

y 

provided g s > Ck/{1 + |A|) S . 
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Proof. Suppose the maximum M is attained at y ^ x. Then 



CkM CM M 



2g s {l + \X\Y ~ 2C 2 ' 
a contradiction. □ 
Corollary 2.2.2. For any s < and g s > Ck/(1 + |A|) S 

(\\Gx{x,x)\\') < ^ 
Proof. By Proposition 12.11 with y = x and Corollary 12.2. 1[ 

(\\Gx(x,x)W) < 2(1 ^ |A|)S {g-K(\\G x (x,x)W) + 1} 
<l(\\Gx(x,x)\\')+ ( 



2 m, 2(1 + |A|) 4 

therefore 

(HG A (x,x)ir>< r 



1 + |A|)' 



□ 



Proof of Theorem [7771 For x = y the inequality follows from Corollary l2.2.2[ For 
x ^ y apply Proposition |2.1l dist(q;. y) times, and then use Corollary 12.2. 1 1 and 
Corollary !2.2.2l to estimate every term. □ 

Proof of Theorem [772] The proof follows that of Theorem ll.il The main modifi- 
cation (apart from replacing || • || with | ■ |) appears in Lemma [2T2l which has to be 
replaced with 

Lemma 2.3. For s < n , ac ! , , there exists C such that 

— 2ka+kq' 

(|G A (m,n)| s |F(n) - X\ s ) > &-\\G x {m, n)| s )(l + . 
The proof is provided at the end of Section [3l 

□ 
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3 Estimates on ratios of polynomials 

Lemma [231 will follow from 

Proposition 3.1. Let fi be a probability measure satisfying the assumptions Al, 
A2. Let a±, ■ ■ ■ , ai, bi, ■ ■ ■ ,b m G C, and let s,r > be such that rm < a and 
q > (si + rm) — - — . Then 

n-=ik-%i* , , nj =1 ( 1+ M) - 

djx{v) ^ 



where the "x " sign means that LHS < C RHS < C LHS, and the numbers 
C, C > may depend on a, q, Cai, C^i, I, m , r, and s, but not on a,j and 6j. 

Proof of Lemma |2"72l First let us show that the statement holds for very small s > 
0; then we shall extend it to all s < i( ^+2q • ^ e sna ^ cons ider the (slightly more 
complicated) case x ^ y. 
For s sufficiently small, 

(\\G x (x,y)\\*U y) 

< (\\G x (x,y)(V(y) - X)\\ s \\(V(y) - X)' 1 ^ (3 .l) 

< (\\G x (x,y)(V(y) - X)\\ 2s ) 1 J^- ) (\\(V(y) - A) _1 || 2s )^ . 
By the Schur-Banachiewicz formula for the inverse of a block matrix^, 



G\(x,x) G x (x,y) 
G x (y,x) G x (y,y) 



V{x) 

V(y) 



-K. 



2kx2k 



where K 2kX 2k is independent of v(x),v(y). Applying the Schur-Banachiewicz 
formula once again, we obtain: 

n ( \ t (\n \ n/f W L kxk (V(y) - M kxk ) Ad > 
G x (x, y) = L kxk (V{y) - M kxk ) =- 



det(V(y) - M, 



kxki 



where L kxk and M kxk are independent of v(y), and Adj denotes the adjugate (= 
cof actor) matrix. 

r ( m/( v n L kxk (V(y) - M kxk )^(V(y) - A) 
G x (x,y)(V(y) - X) = — 



det(V(y) - M, 



kxk) 



2 See Henderson and Searle iTPl for the history of block matrix inversion formulae 
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Therefore every entry of G\(x, y)(V(y) — A) is a ratio of two polynomials Q 1} Q 2 
of degree < k with respect to the variable v(y). For sufficiently small s > 0, 
Proposition 13.11 implies that for any such pair Q%, Q 2 

Hence 



(||Ga(*,I/)(V(i/) - A)!! 28 )^ < kC(\\G x (x,y)(V(y) - X)\\ s ) v{y) . 
Proposition 13.11 also implies that for sufficiently small s 

(\\(V(y)-\yY s )Z<m + \Mr s - (3-2) 
Indeed, using Proposition ^ . 1 1 we first observe that for sufficiently small s 

(W(v(y)A(y) + B(y) + i)(V(y) - A) _1 || 2s )^ < 1.1k , 
Using now the resolvent identity 

(V(y) - A)- 1 = -{i + A)- 1 + (i + X)-\v(y)A(y) + B(y) + i)(V(y) - A)" 1 , 

we establish (13.21) . 

Returning to (13.11 ), we obtain: 

(\\G x (x,y)\\ s ) v{y) < C(\\G x (x,y)(V(y) - X)\\ s ) v(y) (l + \X\)~ S . (3.3) 

To extend this inequality to all s < 4 "^ 2q , we apply Proposition 13.11 once again. 
Every entry in G\(x, y) and G x (x, y) (V(y) — A) is a ratio of two polynomial func- 
tions of v(y) whose degree do not exceed k. By Proposition [3TTJ the expressions 

\Q2{V)\ S 

are comparable as long as q > ^zj^, that is, s < k ^ ka - Therefore (13.31) remains 
valid in this range of s. Averaging over (v(z)) z ^ y , we obtain 

(\\G x (x,y)\n < C(\\G x (x,y)(V(y) - X)\\ S }(1 + \X\)~ S . 

□ 



13 



Proof of Proposition [XT] 

Lower bound. Choose R > so that //[— R, R] > 1/2 (for example, take R = 
max(l,2C A2 y/ q .) Then 



dfi(v) > 



]\k=l\ v - h k\ r J-R 

. „_iri|a j |>2i?( :L + \ a j\) S f R -TT | is t / x 

Now, for any < t < 1, the set |i> — a 3 -| < i} can be covered by / intervals of 
length Ct x l l . Therefore, when t > is sufficiently small, 



| ( \v-aj\ <t > < C 2 t Q/i < 1/4. 



|a,|<2i? 



Then 



/ n \^-^\ s Mv)>t s /^>c^>c^ n 

• R |a J |<2/? |<Xj|<2.R 

Upper bound. Let us start with several reductions. First, it is sufficient to consider 
the case a x = • ■ ■ = aj = a, &i = • • • = b m = b. This follows from the Cauchy- 
Schwarz inequality 

n' iD-fl-i s 1 m r r in n \ si ih 



Second, 



v-b\ rm w ~~ [J |w-6| rm J l^-^l 



so it is sufficient to consider the case a = 0. Third, we can assume that |6| > 1, 
since for \b\ < 1 the regularity condition Al implies 

\v\ sl dfi(v) <c f [ Mv) ! | 6 i 8 f /* d^QQ \ <c < C 6 



w _<j|«n - yj \ v _tyrm-sl ' ' J \ v - b\r™ ) ~ ° ~ (1 _|_ 
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Therefore we need to show that for \b\ > 1 

/\v \ sl dfi(v 
\v - b\ rm 

Let us divide the integral into two parts: 



<C\b[ 



+ 

v\-\b\ \>\b\/2 J|6|/2<|w|<3|6|/2 



By A2, the first integral is at most 

\v\ sl d/i(v) < C 6 \b\- 



2 N rm 



.\b\ 

Let us estimate second integral. 

'3|&|\ s ' f dfi(v 



|6|/2<|b|<3|6|/2 V 2 / y|6|/2<|u|<3]6|/2 \ v - &|' 



< C 7 \b\ sl I /2 \ \b\/2 < \v\ < 3\b\/2 , \v - b\ < r™ \ dt 



oo 





ftfy— sl — rm t*b^ 

\sl 



Cr\b\ sl { + / + 

'0 Jb- sl ~ rm Jb'< 



(3.4) 

where 7 > is a number that we shall choose shortly. The first integral in (13.41) is 
at most t)- sl - rm . The second integral is at most 

6>{M > I&I/2} < C 8 b^- q < C 8 b- Sl ~ rm 

as long as 

7 < q — si — rm . (3.5) 

The third integral is at most 

/i\\v-b\< t~™\ dt<C 9 [ r™ dt < C 10 |6| _7( ™ _1) < C u \b\- sl - rm 
as long as 

sl + rm 

7 > *_ _ 1 ■ (3-6) 

rm 

Since q > (si + rm) J* , we can choose 7 that satisfies both (13.51 ) and (13.61) ; 
then we obtain the claimed estimate. □ 
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Now we prove Lemma [23] 



Proof. We shall prove that 

(|G A (m,n)| s |y(n) - A|% n > ^({G^m, n)\ s ) Bn (l + |A|) S . 

First, 

(|G A (m,n)r/ 2 )| n <(|G A (m,n)nV(n)-A| s ) Bn (|V(n)-Ar s ) Bn) 
and, as above, 

<mn)-Ar> Bn <C(l + |A|r . 
Therefore it remains to show that 

(|G A (m,n)| s ) Bn <C(|G A (m,n)| s / 2 )| n . (3.7) 

For simplicity of notation, let B n = {t>i, ■ • • , Vj} (here 1 < J < k). Cramer's 
rule (or the Schur-Banachiewicz formula) shows that, as a function of every Vj, 
Q(vi, ■ ■ ■ ,vj) = G A (m, n) is a ratio of two polynomials of degree at most k. The 
next multivariate version of Proposition 13 . 1 1 concludes the proof. □ 

Proposition 3.2. Let Q be a function of J variables v i, • • • , v j which, as a func- 
tion of every Vj, is a ratio of two polynomials of degree at most k. Then, for any 
probability measure ji satisfying the assumptions Al, A2 and for any s < k ^" 2 a) ' 



dp,{v i) • • • dfi(vj) \Q(vi, ■■■ ,vj)\ s \ 

2/s 

< Ck, a , q ,c A1 ,c A2 \ I dp,{v x ) ■ ■■dfx(v J )\Q(v 1 , ■ ■ ■ ,vj)\ s/2 



Proof. If J = 1, (13.71 ) follows from Proposition 13. li Then we proceed by indue - 
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tion on J. By case J = 1, 



dn(vi) ■ ■ ■ dfj,(vj)\Q(vi, • ■ ■ ,Vj)\ s 

<d dfi(vx) ■ ■ ■ dnivj-x) \ dfi(vj)\Q(vx,--- ,vj)\ s/2 



C x j dii{v x ) ■ ■ ■ dfi(v 

dp(vj)\Q(v lr • • ,vj)\ s ' 2 I d^v'j)\Q{v ir ■ ■ Xj)\ s/2 



= C\ j dfi(vj)dfi(v'j) J djj,{vi) ■ ■ ■ dfi(v 

\Q(v u --- ,vj)\ s / 2 \Q(v u --- ,v'j)\^ 

By the Cauchy-Schwarz inequality and the induction step, the last expression is 
at most 



C x I dfi(vj)dfi(v'j) 

dn(vx) ■ ■ ■ dfi(vj-x)\Q(vx, ■■■ , vj) | s J dfx(v[) ■ ■ ■ dfi(v'j_x) |C?(v£, • • • , v'j) \ l 
<C 2 J dn(vj)dn{v'j) 

dfx(vx) ■ ■ -dn(vj-i)\Q(v u ■ ■ ■ ,vj)\ s/2 J dfx(v[) ■ ■ • d/i(«j_i)|QK, • • • 

C 2 |y dfi(vx)---d^vj)\Q(vx,--- ,vj)\ s/2 \ . 

□ 



4 Dynamical localisation 

The proof of dynamical localisation is based on the notion of eigenfunction cor- 
relators, introduced by Aizenman in H). 

Let us start with some definitions, which are adjusted from the lecture notes 
of Aizenman and Warzel [3] to our block setting. Let H be an operator acting on 
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? 2 (z d ) <g> C k . For m, n G Z d , the (matrix-valued) spectral measure /i mn is defined 
by 

e CM . 



J <frdn mn = (j>(H)(ia,n) 



The eigenfunction correlator Qi(m, n) corresponding to a finite interval I C R 
(on the energy axis) is defined by 

Q/(m,n) = sup{||0(fr)(m,n)|| | supped, |0| < 1} . 

Obviously, 

sup |e^ 7 (m,n)| < Q/(m,n) (4.1) 

for any t > 0. 

The eigenfunction correlators can be also defined for the restriction of H to 
a finite box A (we denote this restriction by the superscript A). In this case, it 
satisfies the following inequalities (the first one is an equality in the scalar case, 
cf. 0): 

Lemma 4.1. 

Qf(m, n) < foa € - J \\G A +i0 (m, n)\\ l ~*dX < k . 

Proof. For any eigenvalue v of H A , define a k x k matrix 

M v = ip(m) <S> ip{n) '■ u i— > ^^(^(n) ■ M)^(m) , 
where the sum is over all eigenfunctions if) of _?f A corresponding to v. Then 

0(F A )(m,n) = ^0(z/)M,, 

whereas 

V 

(where now the sum is over all eigenvalues of H A .) Therefore 

||0(tf A )(m,n)|| <^||M,|| = Hm| I ||G A (m, n)\\^dX . 
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The equality can be proved by representing / = W I v as a disjoint union of neigh- 
bourhoods of v E I, and noting than 



G^(m,n) = ^- + 0(1), X^v 

V — A 



Also, 



E IWI ^ E IW m ) ® ^(n)|| = E ll^(m)||||V»(n)|| 

< {En^ m )ii 2 E^( n )ii 2 l = k - 

□ 

Now assume that H = H u is a random operator of the form (|1.2I) . where the 
random V(m) are independent and identically distributed, K(m, n) depends only 
on m — n. We shall prove 

Theorem 4.2. Let < s < 1, an J suppose for every box A, every X E I, and 
every m, n G A 

(||G^ (m,n)|| s ) < Cexp(- 7 dist(m,n)) 

/or some C, 7 > 0, where G A is the resolvent of the restriction of H to A. Then, 
for every m, n E Z d , 

(Q/(m,n)) < C'dist M (m,n)exp(-^dist(m,n)) . 

8 

Remark. A similar statement can be proved for potentials of the form rti.iD . 

Proof. We shall prove the estimate in a large box A containing m, n (uniformly 
in the size of A). Let A m and A n be two boxes of radius R = |_dist(m, n)/2j, 
centered at m,n, respectively. According to the resolvent identity, 

G A x+t0 (m,n) = g- 1 E ^(m, w)K(w, w')G A +i0 (W, n) , 

ww'69A m 

where the sum is over all pairs ww' such that w E A m , w' ^ A m , w ~ w'. 
Therefore 

\\G A +iQ (m,n)\\<Cg- 1 max ||G^ (m,w)|| E H G a-boK n)|| . 

ww'69A m 
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Now we apply (8J Prop. 5.1] (which holds in the block-operator setting). It shows 
that, with probability at least 1 — C'R 2d exp(— ^sR/8), one can decompose / = 
-fm U I n so that for every ww' 6 <9A m and A G I m 



max ||G^. (m,w)|| <Cexp(- 7j R/8) 



ww'69A 

and for every ww' 6 <9A n and \ E I n 



'e8A 

Therefore, 



max ||G^ (n,w)|| < C exp(- 7 i?/8) 



lim n - / IIG^o^n)!! 1 -^ < C"^- 1 exp(- J R 7 /8) J R d - 1 



e 

e"ATo2 It 

and the same estimate holds for the integral over I n . Therefore finally 

(Q 7 (m,n)) < C IV g^R*- 1 expi-R-f/8) + C'kR 2d exp{- 1S R/8) 
< C v R 2d exp{- 1S R/8) . 



□ 
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